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Abstract 

We consider the ubiquitous Duffing oscillator driven by a weak periodic force with a time-delay 
linear feedback term. For a random time-delay we numerically explore the effect of strength of 
feedback term (7) and standard deviation a of Gaussian white noise type delay time on resonance. 
We show the occurrence of single and double resonance and bring out the mechanisms associated 
with the resonances. 
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I INTRODUCTION 



Nonlinear oscillators are capable of displaying variety of resonances depending on the 
nature of perturbations. Some of the interesting resonances are nonlinear resonance [1, 2] 
(the perturbation is an external periodic force), Coherence resonance [3, 4] (induced by noise 
in the absence of an external periodic force), stochastic resonance [5, 6] (induced by noise 
in the presence of an external periodic force), ghost resonance [7, 8] (resonance at a missing 
fundamental frequency of the multi- frequency force), vibrational resonance [9, 10] (resonance 
at a low-frequency force induced by a high-frequency force) and chaos resonance [11] (due 
to a chaotic perturbation). In recent years great interest has been focused on time-delayed 
nonlinear dynamical systems and many fascinating dynamics have been observed [12, 13]. 
The time-delay can be a single constant, multiple constants, continuous over a range, state- 
dependent, distributive, periodically varying with time and random. The present paper is 
concerned with stochastic and vibrational resonance with random time-delay. 

In certain conditions the time-delay can be random. For example, in biochemical reac- 
tions, a source of time-delay can be the response time of the last step which can be random 
due to the molecular thermal motion. In population dynamics when the time-delay due to 
the gestation or incubation times depend on the variations of the external environment and 
living conditions due to climate and weather changes then it has to be treated as a random 
one [14, 15]. Random delay enters in the immune response models if the outbreak of a disease 
caused by an infection are due to complicated biological and environmental processes [16]. 
Random delay is found to give rise stabilization of steady state and periodic orbits [17-21], 
synchronization [22, 23], H-infinity filtering [24, 25], state estimation [26], existence of a 
stationary random solution [27], controlling of systems [28] and bifurcations [29] have been 
studied in certain systems with random delays. Random delay induced resonance with and 
without an external periodic force is realized in the Langevin equation [30, 31]. Solutions of 
certain linear differential equations with random delay have been analyzed in ref. [32]. 

The goal of the present paper is to study the occurrence of resonance with random delay- 
time. We choose the reference model system as the under damped Duffing oscillator with 
linear time-delayed feedback and driven by an external weak periodic force. First we show 
that typical stochastic resonance occurs in the system with additive Gaussian white noise 
when the time-delay r is a constant. Next we show the resonance occurs in the absence of 

2 



additive Gaussian white noise but with the time-delay r being Gaussian. We consider the 
effect of 7, the strength of the feedback term and a the standard deviation of the Gaussian 
noise representing the time-delay r. Single and double resonances are found when 7 or a 
is varied. We explain the mechanisms of the resonance using time series plot and phase 
portrait. 



II RESONANCE WITH FIXED TIME-DELAY 



The equation of motion of the ubiquitous Duffing Oscillator with a linear time - delayed 
feedback, driven by a periodic force of frequency uj and additive noise rj is given by 

x + d x + uj^Xi + f3x\ + ^x(t — t) — f cos (uot + rj(t)) (1) 

r)(t) is a Gaussian white noise with mean zero and standard deviation a. r is the time- 
delay and 7 is the strength of the feedback. Throughout our study we fix the values of the 
parameters as d — 0.5, uJq — 1, (3 — 1, / = 0.1 and 7 = 0.1. The potential of the system is a 
double-well form with a local maximum at x = 0 and two local minima at x± = a/| uJq \ /(3. 
Euler method with step size h = 0.001 is used to numerically integrate the eq.(l). Before 
investigating the effect of random time-delay in the absence of r](t), first we consider the 
system (1) with constant time-delay. We fix r = 0.5 and numerically integrate eq.(l). After 
leaving the solution corresponding to first 500 drive cycles as a transient, we calculate the 
response amplitude Q over next 30,000 drive cycles. Q is given by 



Q = VQl + Ql ■ (2a) 



with 



2 r nT 

Q s = — J x(t) sin(27rf/T) dt , (2b) 

2 r T 

Q c = — \ x(t) cos(2vrt/T) dt , (2c) 
nT Jo 

where T = lix joj and = 30000. Figure 1 presents the numerically computed Q as a function 
of the parameter a (characterizing the strength of noise) for two values of uj. This figure 
shows a typical noise-induced resonance curve. For both uj — 0.1 and uj = 1 as the value of 
a increases from a small value Q increases, reaches a maximum value at a critical value of 




FIG. 1: Variation of the response amplitude Q with the standard deviation of the noise r](i) for 
the system (1) for (a) uo = 0.1 and (b) oj = 1. The values of the parameters in the system (1) are 
d = 0.5,^ = -!,/?=!,/ = 0.1, r = 0.5 and 7 = 0.1. 



0 and then slowly decreases. For uo = 0.1, Q is maximum at a — 0.067 while for u — 1, Q 
is maximum at o = 0.021. 

The dynamics of the system is different at resonances for u — 0.1 and 1. In fig.2 the 
evolution of the state variable x is plotted for four values of a for ip = 0.1. In fig. 2a for 
a = 0.001, x(t) is confined to the region x > 0 (there is another trajectory bounded to the 
interval x < 0). x(t) is similar to the input periodic signal, however, it is perturbed by the 
noise term. At a critical value of a noise-induced switching between the regions x < Oand 
x > 0 occurs. This is shown in fig. 2b for a = 0.04. The time intervals between successive 
switchings are not the same. Switching intervals are random. At another critical value of 
er(= 0.067) almost periodic switching between the two wells occurs. This is the case in fig.2c. 
The mean residence time of a trajectory in each of the two potential wells is ~ T/2 where 
T = 2tt/uj. For o » 0.06 rapid switching of the trajectories between the regions x > 0 and 
x < 0 happens (fig. 2d). Figure 3 displays x(t) versus t for uj — 1 and a = 0.021. Though 
the response amplitude Q becomes a maximum at a = 0.021 there is no periodic switching 
of the trajectory between the two wells. This example shows that periodic switching is not 
necessary at stochastic resonance. 



Ill RESONANCE WITH RANDOM TIME-DELAY 



Next, consider the system (1) without the additive noise term but treat the delay time 
r as random. Particularly choose it as a Gaussian white noise with mean zero and variance 
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FIG. 2: x{t) versus t of the system-1 for four values of a with uj = 0.1. 
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FIG. 3: x(t) versus t of the system-1 for a = 0.021 and oo = 1. The dashed line is the rescaled 
input signal / cos ut. 



a 2 . Because a Gaussian random number can be a negative the feedback term is chosen as 
x(t— | r |). The equation of motion of the system is 

x + d x + oj^Xi + (3x1 + ix(t— | r |) = / cos (cut) (3) 

We treat the strength of the feedback described by the parameter 7 as the control parameter. 
Figure 4 presents the dependence of the response amplitude Q on 7 and a when ou — 0.1, 
for a range of fixed values of cr(0 < a < 3) a single resonance when 7 is varied (fig.4a). The 
value of Q at resonance decreases with increase in the value of a. Q at resonance becomes 
relatively much for a > 1.75. As shown in fig.4b a different behavior takes place for u — 1. 
For each fixed value of a there are two resonances when is varied from, say -1. The second 
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resonance is the dominant resonance. 

In fig. 5 we plot jr, the value of 7 at which resonance occurs and Q max , the value of Q at 
7 = 7r, as a function of a for u = 1. The value of 7 at which first resonance occurs (denoted 
by 7^) and the corresponding value of Q (denoted as Qmax) are almost constant with the 

parameter a of the random time delay. 7^ sharply decreases with o and then increases 

(2) (2) 
almost in a linear fashion. In fig. 5b Q max displays sinusoidal type variation. Q ma x increase 

rapidly with a and then attains a saturation. 
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FIG. 5: The variations of (a) 7^ (the value of 7 at which resonance takes place) and (b) Q T 
(7 = 1r) as a function of the parameter a associated with the randomness of the time delay. 
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FIG. 6: x(t) versus t of the system-2 for four values of 7 when u = 0.1 and a = 0.1. The dotted 
curve is the rescaled / cos ut and continuous curve is the numerically computed x(t). 



We understand the mechanism of observed stochastic resonance using trajectory plot. 
Figure 6 shows evolution of x(t) for four fixed values of 7 and for a = 0.1 and u — 0.1. For 
this value of uo the resonance occurs at 7 = 0.6. In all the subplots of fig. 6 the dotted line 
is the rescaled input periodic signal / cos out. In fig.6a (7 = 0.6) far before resonance x(t) is 
bounded in the region x > 0 and is almost periodic with the period of the input periodic 
force. For 7 = 0.6145 (fig.6b) x(t) switches between the positive and negative values, that 
is, jumps between the two wells of the double-well potential. The motion is nonperiodic 
at the time intervals between consecutive jumpings are random. In fig. 6c corresponding 
to the resonance case (7 = 0.68) the solution x(t) appears almost periodic and we clearly 
notice cross-well motion and synchronization between x(t) and / cos ut (the external driving 
force). Compare fig. 6c with fig. 2c when the resonance is induced by the additive Gaussian 
white noise. In both the cases the mean residence time is T/2, however, in the additive 
noise induced resonance case x(t) is not periodic. The erratic switching between the wells 
seen in fig. 2d is not observed for (3 » 0.68 in the random delay. The number of switching 
between the wells per drive period increases with increase in 7 (mean residence time in a 
well decreases with 7), however, the motion is not highly random. 

For uj = 1 as shown in fig.4b there are two resonances when 7 is varied. The two 
resonances occur at 7 = 0.41 and 1.72. Now we account these two resonances. Figure 7a 
displays phase portrait of the system for four values of 7 for a = 0.1 and u — 1. As 7 
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increases from zero the center of the orbit moves towards the origin as the amplitude of the 
orbit increases. For a = 0.41 the amplitude is maximum however the orbit is confined to 
a single well alone. As the value a increased further from 0.41 the amplitude of the orbit 
decreases and the center of the orbit further moves toward the origin. At 7 = 1 the response 
amplitude Q becomes maximum and the center of the orbit becomes the origin. For 7 > 1 the 
center of the orbit remains the same while the amplitude Q shows nonmonotonic variation. 
At 7 = 1.72 Q becomes maximum. This is the second resonance. For 7 > 1.72 the response 
amplitude decreases. The above scenario around the second resonance is presented in fig. 7b. 
The figs 7a and 7b clearly bring out the differences between the mechanisms of random 
delay induced two resonances. For sufficiently small uj, for example uj — 0.1 (fig. 6), the 
observed random time-delayed feedback induced resonance is similar to the typical stochastic 
resonance induced by an additive Gaussian white noise. Here the resonance is due to the 
almost periodic switching of the trajectories between the two wells. For large uj, for example 
uj = 1 (fig. 7) induced double resonance is like the vibrational double - resonance observed 
in underdamped monostable and bistable oscillators induced by a biharmonic force [9,10]. 
There is no periodic switching between the potential wells at resonance. In the random 
delayed system (2) we can introduce an effective potential due to the inclusion of random 
delayed feedback term. Then as introduced in ref. [9,10] we can define resonant frequency 
uj r which depends on the parameters 7 and a in addition to uj^ and (3. As 7 or a varies 
uj r also vary nonmonotonically so that when uj r matches with the frequency uj, Q becomes 
maximum and a resonance results. 

Resonance is found when the parameter a is varied. For uj — 1 and for 7 = 0.4 and 
0.65 the response amplitude Q is computed for a range of values of a and is shown in fig. 8. 
For 7 = 0.4 single resonance at a — 0.68 is observed. When 7 = 0.65 Q is maximum 
at a — 0.2 and 2.65. Resonance at these values of a is also not associated with periodic 
switching between the potential wells. In fig.9 phase portrait of the system (2) is shown 
for three values of a for 7 = 0.4. For all the values of a the center of the orbit is around 
a/| uJqO + 7 I 7/3 ~ 0.775. The orbit corresponding to o — 10 lies in between the orbits of the 
cases 7 = 0.1 and 0.68. That is, the orbit is not continuously expanding. It expands upto 
7 = 0.68 and then shrinks with further increase in 7. We can expect uo r ~ uj at 7 = 0.68. 

In fig. 10 phase portraits for four values of a with 7 = 0.65 are plotted. For a < 0.22 as 
a increases from a small value the center of the orbits remain the same while the size of the 
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FIG. 7: x{t) versus t of the system (2) with to = 1 and a = 0.1 and for four values of 7 around 
(a) first resonance and second resonance. In (a) the values of 7 for the orbits from right to left 
are 0.1, 0.25, 0.41 (resonance case) and 1 (at which Q is maximum). In (b) the values of 7 for the 
orbits 1, —4 are 1, 1.7, 1.72 (resonance case) and 2 respectively. 
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FIG. 8: Variation of the response amplitude Q with the standard deviation of the noise rj(t) for 
the system (2) for two values of 7 = 0.1 and 7 = 0.65. The values of the parameters in system (2) 
are d = 0.5, = -1, 0 = 1, / = 0.1, and r = 0.5. 



orbit ( and the value of Q ) increases and becomes maximum at a = 0.2 ( first resonance 
value) and then decreases. At a = 1.22 onset of cross- well motion occurs. This is shown in 
fig. 10b for a = 1.25. At a — 2.65 Q becomes a maximum (second resonance). In fig. 10 we 
can clearly notice difference between the orbits shown in fig. 10c and those shown in figs 10b 
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FIG. 9: The phase portrait of the system (2) for three fixed values of cr. Here u> = 1 and 7 = 0.4. 
Q is maximum at a = 0.68. 
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FIG. 10: x versus x of the system (2) for u = 1, 7 = 0.65 and for four fixed values of cr. 



and d. The orbit of the system is a cross-well orbit in figs 10b - d, however, Q is maximum 
for cr = 2.65. Therefore we can say that u r = u at 2.65. 



IV CONCLUSION 



In the present work we have reported the resonance induced by the random time-delayed 
feedback in Duffing Oscillator. For small values of u, the observed resonance is like an ad- 
ditive noise induced stochastic resonance. That is, single resonance occurs and at resonance 
almost periodic switching between the two wells takes place. For large values of u double 
resonance is also realized. At the first resonance cross-well motion is not found to occur. 
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On the other hand, second resonance takes place far after the cross-well motion. 

As the random delay induced resonance in the presence of external periodic force is 
not strictly like a stochastic resonance. It is important to investigate the occurrence of 
other types of resonances like coherence, vibrational and ghost resonances in the presence 
of random delay. Further the study of various resonances with distributive delay, state- 
dependent delay and time-varying delay can bring out interesting and novel dynamics. 
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